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REVISION NOTICE 

This description replaces previous descriptions of "Matrix Multiply 4," 
program Dl-239. 2. Program references have been changed to current 
designations. 



FUNCTION 

"Matrix Multiply 4" enables the interpretive system to multiply 2 
matrices A(i x j) and B(i x j), where ja and ib are the same but the 
matrices are not necessarily square. The routine is entered and left 
in machine language, but uses the Floating Point Interpretive System 4 , 
program Hl-24.3 for all arithmetic. 



INPUT 

The elements of matrix A begin in Ao and the elements of matrix B 
begin in Bo. The data are stored consecutively in double precision 
floating point format. 



OUTPUT 

The elements of the product matrix C are stored in double precision 
floating point format beginning in Co. 



SC 0070 



1 of 2 



PRINTED IN U.S.A. 



MATRIX MULTIPLY 4 



CALLING SEQUENCE 

Location Order 

a R 

a + 1 U 

a + 2 Z 

a + 3 Z 

a + 4 Z 

a + 5 Z 

a + 6 Z 

a + 7 Z 



Address No t e s 

Lo Initial location of "Matrix 

Multiply 2." 
Lo Initial location of "Matrix 

Multiply 2." 
Io Initial location of interpretive 

routine, 
iaja ia in track, ja in sector, 
ibjb ib in track, jb in sector. 
Ao Location of matrix A. 
Bo Location of matrix B. 
Co Location of product matrix C. 



TIME 



Approximately 4(ia)( ja)( jb) milliseconds are required 



STORAGE 

2 tracks, 32 sectors (160 words) of memory are required for this sub- 
routine plus storage for matrices A, B, and C. 



PROGRAM STOPS 



Location 



Lo + 35 



Meaning 

ja not equal to ib. Do not continue. 
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Job #OlU9 



DOUBLE PRECISION FLOATING POINT MATRIX MULTIPLICATION 



FUNCTION : 

To multiply two matrices A(i x j) and B(l x j) where ja and ib are the same but 
the matrices are not necessarily square. The routine is entered and left in machine 
language, but uses DPFP for all arithmetic. 

INPUT : 

The elements of matrix A beginning in Ao and the elements of matrix B beginning 
in Bo. Data are in DPFP format. 

OUTPUT ; 

The elenents of the product matrix C in DPFP format stored consecutively 

beginning in Co. 



CALLING SEQUENCE: 



Location 


Order 


Address 


a 

a + 1 
a + 2 
a ♦ 3 
a + U 

a + 5 
a + 6 
a + 7 


R 
U 
Z 
Z 
Z 
Z 
Z 
Z 


Lo 

Lo 

Lo of DPFP 

iaja 

ibjb 

Ao 

Bo 

Co 


TIME: 






Approximately k iajajb 


seconds. 


STORAGE: 







Notes 

■ m il ■ ii | ■ i n 

Initial location of this subroutine 
ii tt ii it n 

Initial loe, of DPFP 
ia in track, ja in sector 
ib in track, Jb in sector 
Location matrix A 
Location matrix B 
Location product matrix C 



2 l/2 tracks for this subroutine plus storage fop matrices A, B, and C. 
PROGRAM STOPS : 

Location Meaning 

Lo + 3£ ja not equal to ib. Do not continue, 
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